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ABSTRACT 

The structure of molecular cloud cores supported by thermal pressure and a poloidal magnetic 
field is reinvestigated in the magnetohydrostatic and axisymmetric approximation. In addition to 
oblate configurations found in earlier work, solutions yielding prolate spheroidal shapes have also 
been obtained for a reference state described by a uniform sphere threaded by a uniform background 
magnetic field. The solutions for prolate configurations are found to be relevant for lower masses 
than for their oblate counterparts. Of particular importance is the result that the prolate cloud cores 

/ \ 2 / \ -1/2 

have radii less than a maximum given by 0.25pc ( 2km/ s ) ( lo-^dyne/cm'* ) > wnere a is the sound 
speed and P ex t is the external pressure of the background medium. The existence of such solutions 
obviates the presence of toroidal fields in such modeled structures. 

Subject headings: star formation, molecular cloud cores, magnetic fields 



Oh- 
i < 

O 



> 

O 

m 

o 



1. INTRODUCTION 

The study of dense molecular cloud cores are of 
paramount in terest because they are the sites where 
stars form fsee lShu et. ldll987l ; lMcKee fe Ostriker!l2007l . 
for reviews). In situations where other physical ef- 
fects in comparison with gravity and thermal pressure 
can be neglected, these cores have been modeled as 
isothermal spheres in hydrostatic equilibrium bounded 
by an ex t ernal pressure, kno wn as Bonnor-Ebcrt spheres 
(|Bonnorl [1955 lEber^ [l955h . For some objects, such 
a model offers remarkable agreement with observation 
(lAlves et. al.ll200HlEvans et. al.ll2001l; 



iSchnee fe Goodmanll2005l ; iStutz et. al 



Kirk et. al.|[2005l ; 
[2007D However, 



surveys of large samples of dense cores in dark clouds 
reveal that spherically symmetric cores are the excep- 
tions rather than the rule. The projected aspect ra- 
tios can significantly differ from unity, making the in- 
trinsi c geometry even m ore elongated ([Myers et. al. I 
[19911: IJiiina et. all H999h . If we assume these cores 
are axisymmetric and are randomly oriented in the 
sky, statistical analy sis su ggests many of them are pro- 
late (|Mvers et. al. Ml991t iRvdenl Il996f) . On the other 
hand, oblate shapes may fit better with the observed 
distribution of shapes for intrinsically tr iaxial cores 
(|Jones. Basu. fc Dubinskl 120011 ; lTassisll2007[ ). 

The formation of prolate cores presents certain theoret- 
ical challenges. If the formation history of these cores is 
dominated by quasi-static contraction regulated by am- 
bipolar diffusion, t he resulting geometr y is a sequence 
of oblate sphe roids (|Mouschoviaslll976aHbl ; iNakanoll 19791 ; 
ILizano fc Shu1ll98l . This is a natural outcome since 
matter can contract freely along the field lines, but 
must overcome the additional magnetic pressure and 
tension perpendicular to them. The inclusion of rota- 
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tion produces even flatter morphologies, but there is lit- 
tle observational evidence supp orting the view that ro - 
tation is i mport ant for cores (I Goodman et. al.lll993f ). 
iTomisakal (fl99l[ ) and iFiege fc Pudritzl (|2000f ) invoked a 
helical magnetic field to explain the prolate cores and 
were able to produce a large range of aspect ratios con- 
sistent with observation. Here, the hoop stress generated 
by the toroidal field component facilitates the confine- 
ment of matter toward the axis. However, the presence of 
toroidal fields requires a current flowing along the sym- 
metry axis, which is diffi cult for astrophys i cal sy stems 
to realize. We note that iCurrv fc Stahlerl (|2001l ) were 
able to obtain both prolate and oblate cores with only 
a poloidal field by specifying the shape of the core and 
solving for the mass-to-flux ratio after an equilibrium 
had been found. In this case, the mass-to-flux ratios are 
distinct for the two types of cores. 

In this Letter, we reexamine the structure of molecu- 
lar cloud cores, which are treated as magneto-hydrostatic 
equilibria, with the ultimate goal of providing a grid 
of models that can be used for interpreting prc-stellar 
core density structures. We have found that magnetized 
structures for both prolate and oblate configurations can 
be produced for a wide range of aspect ratios from a 
purely poloidal field configuration with a generic func- 
tional form of the mass-to-flux ratio. This is in contrast 
to previous work where prolate configurations were pro- 
duced only with a toroidal field or a specific mass-to- 
flux ratio. The findings of this work help to elucidate 
the physical parameters that dictate the morphologies 
of dense molecular cores. The theoretical framework in- 
cluding the assumptions and equations of the model are 
described in §2. The method of solution and the numer- 
ical results are presented in §3. Finally, we summarize 
and discuss the results in §4. 

2. FORMULATION AND BASIC EQUATIONS 

For mathematical tractability, the figures of equilib- 
rium are restricted to be axisymmetric and described in 
cylindrical coordinates (tu, z). Consider a cloud core 
described by an isothermal equation of state P = a 2 p, 
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where P, p, and a are the pressure, density, and a con- 
stant sound speed. Let this cloud be embedded in an 
environment characterized by a large scale background 
magnetic field Bo = B$i and an external pressure P e xt- 
Pressure equilibrium requires the cloud to have density 
po = P cy ±/cL 2 on its surface. From these parameters, 
along with the gravitational constant, G, the units of 
length, density, pressure, gravitational potential, mass, 
magnetic field, and magnetic flux are taken as 



l ° - yiwW' 

M = An poll 



po, 
-Bo, 



P, 



ext i 11 ) 

$o = 2irB ll 



respectively. For typical cores, if we take a = 0.2 km 
s _1 , P C xt = 10" 12 dyne cm -2 , and B = 20 pG, then 
l = O.lApc, po = 1.5 x l0 3 m H cm" 3 , M = 1.3 M Q , and 
$o = 2.5/iG pc 2 , where ran is the mass of a hydrogen 
atom. In these units, the magneto-hydrostatic equilibria 
are governed by Poisson's equation for the gravitational 
potential, V, as 

V 2 V = o, 



the divergence free condition of the magnetic field 



VB = 0, 



1) 



(2) 



and the force equation 







-Vlogp — W + (/3p)" 1 (V x B) x B. (3) 



In the last equation, (3 = Anp^d 2 /B 2 . The other two 
equations of magnetohydrodynamics (MHD), the conser- 
vation equation, and the induction equation are identi- 
cally satisfied for d t = v = 0. 

The divergence free condition on the magnetic field 
can be satisfied using the flux function <!> = J B • e z ujdw 
(recall a factor of 2tt has been absorbed into the unit 
of flux) . An axisymmetric and purely poloidal magnetic 
field can be uniquely specified by the flux function as 



B = w x e„ 



(4) 



In this case, B is everywhere tangent to the contours 
of <&. In terms of the flux function, the force equation 
becomes 

= - V log p - W - (/3//T 1 V • (vj- 2 V$) V$. (5) 

The force equation can be projected along the magnetic 
field to obtain B • V/i = 0. where 



h = logp + y, 



(6) 



is the specific enthalpy. This implies that h(vj, z) = h(&) 
is a function of <f> alone. The component of the force 
equation perpendicular to the magnetic field, known as 
the Grad-Shafranov equation, governs the spatial distri- 
bution of field lines, and it can now be manipulated to 
read 

V • (uj- 2 V<f) = -Pp^. (7) 

To close the system of equations, we impose the integral 
constraints that the mass in each flux tube is conserved, 
and is given by 



f z ^ dm 
Jo 



dz 



dm 



(8) 



In the above equation, the z integral is performed over 
constant and dm/d<& is the known differential mass- 
to-flux ratio, obtainable from either observation or an 
evolutionary calculation. The core surface, described by 
Z((f>), is a free internal surface of the problem. To deter- 
mine the location of the core boundary and the specific 



because 



1 on 



z=Z 



enthalpy, we note that h = V 

the surface. Once V is known, both h and Z can be 
determined by solving the equation 



z=Z 



dm 



e w-^—dz 
9$ 



(9) 



along each flux tube. 



3. NUMERICAL METHODS AND RESULTS 

Our numerical scheme for constructing solutions is 
an iter ative procedure similar to th at of IMouschoviasI 
(|1976aD and iTomisaka et. all (|1988aD . For definiteness, 
we adopt a differential mass-to-flux ratio correspond- 
ing to a reference state consisting of a uniform sphere 
threaded by the uniform background field. If the total 
mass of the core is M c , and total trapped flux is $ c , then 



dm 



3A/ e 
2$ c 



1 



(10) 



This configuration was termed "parent cloud" by the au- 
thors, and serves as the initial trial solution. Starting 
with the parent cloud, h and Z(<S>) are computed from 
equation (J9]), which allows a calculation of the source 
functions in equations |T]) and 0. These equations are 
then solved by successive over relaxation to obtain inter- 
mediate solutions. A corrected solution is constructed 
by under relaxation and updated with the next itera- 
tion. The procedure continues until the fractional error 
is reduced to a pre-specified level. Once a converged so- 
lution is obtained, higher accuracies are then achieved by 
a series of mesh refinements. 

We were able to re pro duce the solutions presen ted by 
IMouschoviasI (|1976bD and ITomisaka et. all (|1988bf ) using 
their input parameters. We compare our result visually 
to theirs and estimate the fractional error to be within 
5%. In particular, for a fixed total trapped flux, there 
is a maximum mass beyond which no equilibrium solu- 
tion exists. The numerical yalue of this max i mum m ass 
agrees with earlier work by ITomisaka et. alj (|1988bf ) to 
within 3%. This mass is equivalent to the Bonnor-Ebert 
mass, the maximum stable mass for a given temperature 
and external pressure, but modified by the presence of a 
magnetic field. 

In the magnetically regulated star formation paradigm, 
a key parameter is the total mass-to-flux ratio, A = 
/3 1/2 M C /$ C (or in conventional units A = 2ttG 1/2 M c /$ c ). 
The supercritical clouds (with A > 1) are capable of 
continued contraction leading to star formation, while 
the subcritical clouds (with A < 1) are not. Through 
a process of natural selection, modern day molecular 
clouds are most li kely to be in a marginally critical state 
(|Shu et. al.ll2004D . Clouds with A > 1 would have col- 
lapsed to form stars, while those with A <C 1 would have 
evolved to the diffuse interstellar medium. These the- 
oretical arguments are consistent with the findings of 
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iTroland fc Crutcherl (|2008f) who estimated A - 2 in the 
core s. Thus, we argue that inst ead of fixing the total flux , 
as in lMouschov ias (1976bT) and lTomisaka et. "all ( 



a more convenient parameterization would be to vary the 
total mass while maintaining an order unity value of A. 

The solution space is three dimensional, parameter- 
ized by A, (3, and M c . However, for values of A £ [0.1, 7] 
and M c < 20, the solution is approximately degenerate. 
With a = A//3 1 / 2 = M c /$ c held fixed, varying (3 by 
two orders of magnitude only introduces less than a 1% 
change to the equilibrium solutions, indicating the rel- 
ative insensitivity of the magnetic field strength in this 
parameterization. Wc shall thus use a and M c as our 
basic parameters when discussing the numerical results. 

Of primary interest in this study is the core shape, and 
the aspect ratios of the cores are illustrated in Fig. [T]as 
functions of M c for several values of a. Here, the aspect 
ratio is defined as the ratio of the radial extent to the 
vertical extent of the cloud core surface. An important 
feature to note is that for fixed a, there exists a critical 
mass -/Vf cr it below which the core takes a prolate shape. It 
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Fig. 1. — Core aspect ratios as functions of cloud mass for A = 1, 
and a = 0.5, 0.75, 1.0, 2.0, 3.0. The dashed line indicates the critical 
mass at which the aspect ratio is unity. 

can be seen that the prolate configurations are restricted 
to low core masses. Specifically, for a = 1 and typical 
core parameters (see above), the critical mass is 0.96 M Q . 
Furthermore, as A increases for fixed f3 or as (3 decreases 
with fixed A, the critical core mass increases. The density 
and magnetic field structure for a typical prolate solution 
is displayed in Fig. O 

The solutions for prolate cores are also found to exhibit 
a maximum radius, beyond which such solutions do not 
exist. In particular, the parent cloud radius, R c and 
the mass M cr i t that leads to a final core of unit aspect 
ratio are shown as functions of a in Fig. [3] Here, the 
final core radius is virtually indistinguishable from that 
of the parent cloud. It is evident that the radius is a 
nonmonotonic function of a, in contrast to the core mass, 
exhibiting a maximum at a radius of 1.75 (or about 0.25 
pc for typical core parameters) for a ~ 3. 

4. SUMMARY AND DISCUSSION 
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Fig. 2. — A typical prolate core solution is illustrated in cylindri- 
cal coordinates (to, z). The solid curves are density contours, the 
dashed curves represent magnetic field lines, and the dash-dotted 
curves are contours of the gravitational potential. For this partic- 
ular solution, A = 1, a = 2, and M c = 1. It has an aspect ratio of 
0.582, a central density of p c = 2.5, and a radius of R = 1. 



We have reexamined the structure of dense molecular 
cloud cores supported, in part, by magnetic forces and 
thermal pressure. In most of the parameter space that 
we have surveyed, our results agree with previous studies 
that exist in the literature. Of particular interest is the 
discovery of solutions for which the cores assume pro- 
late shapes. As we have shown (see Fig. [T]), for a fixed 
mass-to-flux ratio, prolate configurations are more rel- 
evant for smaller masses, while the oblate ones are for 
larger masses. 

The formation of prolate cores can be understood in an 
evolutionary sense not too different from that of oblate 
ones. For our choice of the reference state, and a fixed 
a, the core mass is related to the density, pi, by 



ay f 3 



(11) 



Cores with sufficiently small mass represent over dense 
regions, and they must expand to achieve equilibrium. 
Because of the additional external pressure exerted by 
the magnetic field, the gas experiences less resistance in 
the direction parallel to the field lines than perpendicu- 
lar to them. The subsequent evolution naturally leads 
to a prolate figure. Instead of an hour glass shaped 
magnetic field commonly seen acco mpanying oblate cores 
in either theoret i cally c alculations (iMouschoviasI Il976bl; 
i Tomisaka et. al.1 Il988bl ) or observations |Girart et all 
l2006Ll2009n . the magnetic field lines bow outwards in the 
midplanc. However, for our particular choice of the par- 
ent cloud, because the overall magnetic energy is much 
larger than thermal energy for small masses, the outward 
bowing of the field lines is not pronounced in general. Of 
all the prolate shaped models we have constructed, the 
fractional distortion of the magnetic field lines only oc- 
curs at a 10~ 4 level. Other forms of dm/d$ may lead to 
more drastic o utward bowing of the fi eld lines for prolate 
cores (see e.g., ICurrv fc S tabler 200lt). 

It is instruc t ive to i nquir e why previous studies, su ch as 
IMouschoviasI (|1976bt ) and lTomisaka et. all (|1988bt ). did 
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Fig. 3. — Parent cloud radius (solid curve) and core mass (dashed 
curve) for solution with unit aspect ratios. 

not find such kind of prolate cores with just poloidal 
fields, since they have surveyed a large portion of the 
parameter space spanning many orders of magnitude in 
both (3 and central density p r . The answer lie s in the 
total flux trapped by the core. iTomisaka et. al.l (|1988bl ) 
used a uniform sphere with various mass-to-flux distribu- 
tions as a parent cloud to initiate their iterative calcula- 
tion. Our particular choice of parent cloud is threaded by 
a uniform field, which corresponds to their A = 1 model. 
Its radius is related to the trapped flux by R c = (2$ c ) 1 / 2 
in dimensionless variables. As we have shown in §3, all 
prolate cores must evolve from parent clouds of radius 
Rc < 1.75, but other authors chose radii larger than 2.4. 
Because mass and flux scale differently with radius, a 
sufficiently large parent cloud with fixed a will always 
be dominated by gravity, and hence evolve to an oblate 
configuration. Whether adopting a different mass-to-flux 
distribution within the cloud will change this conclusion 
remains to be examined. 

The fact that there exists a maximum radius, i? m ax, for 
prolate cores has interesting observational consequences. 
In conventional units, this critical radius is given by 



-Rmax = 0.25pc 



Po 



0.2km/sJ \10 12 dyne/am? 



-1/2 



(12) 



For typical molecular c loud cores, this value is co nsistent 
with t hose studied bv iMvers et. al. I (|1991h and iRvdenl 
(fl99fih . We may infer more detailed parameters of the 
cores if the properties of the environment can be con- 
strained by observations. In particular, if we are able to 
obtain both the values of A and 



(3 = 0.13 



^ext / Bp V 

10" 12 'dyne /cm 2 \IQ fiG J 



then we may compute the value of a. Inspection of Fig. 
[3] provides an upper limit on the mass of a prolate core. 
As remarked earlier, the masses for prolate configurations 
are predicted to be lower than for oblate configurations 
and, hence, such cloud cores are expected to be charac- 
terized by a lower visual extinction. 

Our analysis is restricted to axisymmetric models with 
only poloidal fields. These assumptions are justified if the 
cores were formed via gradual con traction regulated by 
ambipolar diffusion (however, see iBasu fc Ciolekl 120041 ; 
ICiolek fe Basul 120061 for discussions on the growth of 
non-axisymmctric modes in sheet-like structures). Re- 
laxation of these restrictions yields additional possibili- 
ties, which may shed light on the initial conditions and 
the formation process of the molecular cores. In order 
for magnetic confinement to yield prolate cores, the par- 
ent cloud must be over dense and try to expand ini- 
tially. How to achieve such an initial state is an in- 
triguing question. For instance, numerical simulations 
suggest that collisions between turbulent clouds would 
leave behind triaxial dense cores that tend to be prolate 
(iGammie et al.l2003t iLi et al. 1120041 ; lOffner fe Krumhold 
l2009f h One could consider the possibility that the col- 
liding clouds created the high density region in the first 
place. After the turbulence had decayed, the subsequent 
expansion against an anisotropic external magnetic pres- 
sure may facilitate the formation of the low mass prolate 
cores. We reserve such a consideration for a future study. 
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